Relaxation time ansatz and bulk viscosity of hadron matter 
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The bulk viscosity is calculated for hadron matter produced in heavy-ion collisions, being described 
in the relaxation-time approximation within the relativistic mean-field-based model with scaled 
hadron masses and couplings. We show how different approximations used in the literature affect 
the result. Numerical evaluations of the bulk viscosity with three considered ansatze deviate not 
much from each other confirming earlier results. 
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Introduction. Recently, interest in the transport co- 
efficient issue for hadronic and quark matter has essen- 
tially increased due to clarifying the role of viscosity in 
extraction of flow parameters from heavy-ion collisions, 
see review-article Viscosity coefficients in a weakly 
coupled scalar field theory at an arbitrary temperature 
can be evaluated directly from the first principles using 
expansion of the Kubo formulas in terms of ladder di- 
agrams in the imaginary time formalism 0]. Unfortu- 
nately, similar analysis for more general cases is unavail- 
able. Therefore, in order to evaluate transport coeffi- 
cients in multi-component systems with a strong coupling 
between species, one often uses a kinetic approach. Thus, 
one can exploit the relaxation time approximation to 
the Boltzmann-like quasiparticlc kinetic equations. The 
shear and bulk viscosities of the hadron and the quark- 
gluon plasma phases of strongly interacting matter at fi- 
nite temperature and baryon density were evaluated, see 
Refs. [JQ. In 0-11], a similar analysis was performed 
for purely gluon matter. Besides the use of the relax- 
ation time approximation, one needs to do some extra 
assumptions in order to proceed further. In Q, we stud- 
ied how different ansatze used in the literature affect the 
result for shear and bulk viscosities for gluon matter. Wc 
found that the result for the shear viscosity was robust 
with respect to different ansatz reductions, whereas the 
value of the bulk viscosity significantly depends on them. 
In this note we study how different approximations used 
in the literature affect the result for the bulk viscosity of 
the hadron matter. 

Model equations. As in we describe the hadron 
phase in terms of the quasiparticle relativistic mean- 
field (RMF) based model with the scaling hadron mass- 
couplings (SHMC) successfully applied earlier to the de- 
scription of heavy-ion collision reactions 0, which is 
an extension of the model applied there for the cold 
dense hadron matter. Bearing in mind an application 
to heavy ion collisions we deal with a RMF-based model 
of iso-symmetric non-equilibrium hadron matter with <j- 
and w-meson mean fields and in contrast with ordinary 
RMF models we assume that not only baryon but also 
other hadron masses might depend on the er-mcson mean 
field. Also excitations emerging from the a and wq mean 



fields are incorporated. We study the system with zero 
net strangeness and we use the same hadron set, as in 
[3, • Considering small deviations from local equilib- 
rium we keep only first-order gradient terms. Further 
details can be found in mentioned works [J, Q . 

We start with general expressions for the 
baryon/strangencss 4-currcnt and the energy- momentum 
tensor densities 
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ta > ta are t ne baryon and strange charges of the a- hadron 
(antiparticles are included), d a is a degeneracy factor, 
to*(ct) is an effective mass, f a is a quasiparticle distribu- 
tion function, a and ujq are mean scalar and vector meson 
fields, u" = (w ,0), A° ~w , 
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^nl(c) is the non-linear potential of the a field. 

Conservation laws of the baryon/strangeness 4-currcnt 
and of the energy-momentum tensor densities are read as 
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We assume that the quasiparticlc distribution func- 
tions f a obey a set of kinetic equations 
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where St/ a is the collision term for the given species sat- 
isfying the conditions 



dr a st/ a = o, 
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dr a e a St/ Q = 0, (7) 
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with the quasiparticle energy e a = E a + X®. 

The collision term is zero for the local equilibrium dis- 
tribution, St/^ cq - = 0, where 

= | e K- 8q -P - -«(^)-A':(^)]/T(x'') ±1 |" 1 ) ^ 

with + for fermions, — for bosons, 

E^ = E a (a he *) , 
^(i") = t^ B (x^+t s af ,s(x^-X° a (a^;J ^) , 

CT l.eq. = c7 l.eq.( T)/iJ3;/is ) ) ^ = ° q ' (T, (1 B , (J-s), 

P"B i Ms are the baryon and strangeness chemical poten- 
tials, and the four- velocity of the frame is it M ~ [1, u(x^)] 
for |u| < 1. 

Applying the second Eq. ([5]) for v = and ([7]), and 
using (|2|) we derive the self-consistency conditions: 
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For the equilibrium system these equations coincide with 



the conditions of maximum pressure Sr— = 0, £fL = fj, 

1 CftUo ' OCT ' 

where pressure P = ^T^' cq ' . Latin index i = 1,2, 3. 

In the general case it is impossible to solve the Boltz- 
mann kinetic equations for the strongly interacting multi- 
hadron system appearing in the course of heavy ion col- 
lisions. However the collision term is greatly simplified 
in the so called relaxation time approximation, or better 
to say, in the relaxation time ansatz. We will use the 
expansion near the local equilibrium state 

Stf a = -Sf a /r a , 5f a = f a -f l a ^-(E l a ^-), (10) 

where r a are in general the energy dependent quantities, 
i.e. r a — T a (E l a cq -). These values can be evaluated from 
the cross sections of particle-particle interactions. 

Since the expression for the shear viscosity, rj, is easily 
recovered [|[ within the relaxation time approximation 
and one needs no extra assumptions for that, below we 
study only the bulk viscosity, 

Bulk viscosity. The bulk viscosity is defined as the 
coefficient entering into the variation of T" in the local 
rest frame: 
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and variations of the baryon/strange charge and the en- 
ergy density should satisfy the so called Landau-Lifshitz 
conditions 

sj b,s = E^' S / <w a 6f a = o, (12) 

ST oo = J dT a e a 6f a = 0. (13) 



From the Boltzmann equations within the relaxation 
time approximation we find 
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and e = T^ q \ cf. @, 

Ansatze of Refs. [3|, 0]. Following the line sketched 
in Ref. 0], performing variations in Ref. [|[ we did 

not vary quantities which may depend on the distribu- 
tion function only implicitly, such as E a . This approx- 
imation is well satisfied for non-relativistic systems, see 
Q . Although the validity of this approximation becomes 
questionable in the application to relativistic systems, 
its use allows one to essentially simplify calculations for 
the bulk viscosity, which is important in the case of a 
complicated system of many strongly interacting particle 
species. Therefore, in [3] we considered it as an additional 
ansatz. Then the expression for the value ST" looks very 
simple 
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since E a values are not varied. Using (|16[) we arrived at 
the following expression for the bulk viscosity: 
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Here we introduced the notation 
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As in 0] , we have just assumed the validity of the Landau- 
Lifshitz conditions ilty ) and \1S\) . Using the latter condi- 
tion we found in Q the expression for the bulk viscosity 



CRef.Q 



— Y 

3T ^ 



TaQaiP'' 



,* - 



E a 



X 



.(19) 



One should note that relations (fT2|) . (TT^J) might not be 
fulfilled, until some additional conditions were not im- 
posed, see below. Since we did not impose these extra 
conditions, the use of the Landau-Lifshitz conditions can 
be considered as an additional ansatz. 
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Ansatz of the given work. We again consider the 
relaxation time ansatz. But now we avoid additional two 
assumptions used in [U, H[ . To derive the expression for 
the bulk viscosity, we follow the procedure sketched in 
for gluons (see Sec. III. A there). 

First, from ([9]), (fT0| we find variations of mean fields 
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If these relations are not fulfilled with the particular dis- 
tribution (fT4")h we may still fulfill them doing the shift 



TaQa(p 2 ) -> r a Q a (p 2 ) + y B t B a + y s t s a 



X Sf, 



(27) 



where x and y B ' are some constants. These constants 
are associated with the conservation of the energy and 
the baryon and strange charges. Values of y B and y are 
similar to baryon and strange chemical potentials. Since 
Us 7^ even for hadron matter with zero net strangeness, 
we cannot exclude the term y s . 

If one considers only elastic scattering of particles de- 
scribed by the exact Boltzmann collision term, the re- 
placement (|2T|) generates new solutions of the original 
Boltzmann equation, see [f|. However, one can show 
that for multi-particle systems considered within the re- 
laxation time approximation even with energy-averaged 
values of r a , the above replacement does not result in 
new solutions. Even for one species but with the energy 
dependent relaxation time r(i? l eq -) the replacement does 
not generate new solutions. Thus, we actually fulfill the 
Landau- Lifshitz conditions at the price that the solutions 
of the Boltzmann equations with the collision terms 
are spoiled. The new collision terms are 
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n a = J dr a / Q is the number density of the particle 

species "a" and p s a ° = J dT a ^ f a is the scalar density. 
All integrals in matrix A are calculated in the rest refer- 
ence frame of the fluid with the local equilibrium distri- 
bution function. With the help of these expressions the 
variation ST" can be expressed as 
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In the relaxation time approximation, which wc exploit 
in the given paper, using (|14[) we present the Landau- 
Lifshitz conditions (|T2"j) . (fTB")) as 
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Therefore, the requirement of fulfillment of the conditions 
(|25l) . (|26l) should be considered as an ansatz. 

After performing the replacement (|27p in the condi- 
tions (|25jl , (|26|) , we arrive at the system of linear equa- 
tions for x and y. Finally, we obtain 
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Ansatz of Ref. Q. Above we assumed that the 
original Boltzmann collision terms St/^' eq "(i£*" eq ') = 0. 
However, at calculating the viscosity coefficients in the 
quasiparticle Fermi liquid theory one often uses that 
the original Boltzmann collision terms Stf a ' eq (E a ) = 
also for E a , being functionals of exact non-equilibrium 
distribution functions. This is so because the energy con- 
servation (5-function pre-factors depend on exact particle 
energies. Thus, in the relaxation time approximation 
one can write 

St/„ = -6fa/T a , Sf a = f a -f^{E a ), (31) 

and thereby 



Sfa = Sf a + ^ 

oe a 



de a 91" 
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(32) 



The relaxation time T a is in general different from r a 
introduced above. Note that due to the smallness of 5f a , 
one can consider here r a = T a (E 1 a cq -) as a function of the 
energy E a oq -. 

The difference between the approach Q and that ex- 
ploited in section " Ansatz of the given work" is that fol- 
lowing Q we now express all variations through Sf a and 
r Q rather than through 5f a and r a . Since E a are now 
fixed and thus not varied, we get the same expression 
for the value 5T 11 as ([T5]) but now with the substitu- 
tion 5 f a instead of 8f a . Here one should note that in 
Section "Ansatze of [3,111" the quantity E a [f a ] was not 
varied according to our ansatz, while now the expression 
for ST" [Sfa] becomes a fully correct relation. Within the 
relaxation time approximation using (|31[) we arrive at ex- 
pression p7|) with the only difference that r Q should be 
now replaced by f a . Thus, performing similar calcula- 
tions to those we have done in the previous section we 
rewrite now the Landau-Lifshitz conditions as 
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and solving the corresponding system of linear equations 
for x, y B and y s we find 



CchK = 7p ^2(TaQl(p)) 



(36) 



We note that, as in section "Ansatz of the given work", 
we actually fulfill the Landau-Lifshitz conditions at the 
price of spoiling the solutions of the Boltzmann equations 
with the collision terms £571). 
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FIG. 1: (Color online) The ratio of the bulk viscosity to the 
entropy density as a function of temperature at /is = 0. 



Eq. ([3"6"| formally coincides with that obtained in Q 
but new terms depending on fiB,s ar e involved, cf. (|15[) . 

Numerical results. Details of calculations of the re- 
laxation time T a and of shear and bulk viscosities in our 
SHMC model can be found in Q . Since in such a compli- 
cated multi-particle system, as we study, values r a cannot 
be calculated microscopically but can only be evaluated 
using empirical values of the cross sections, we cannot 
distinguish r a and r a and therefore, as in Q, we consider 
them to be the same. 

In Fig. [JJ we show the bulk viscosity to the entropy 
density ratio at /is — as a function of the temperature. 
The solid line is the calculation of the given work per- 
formed following Eq. (|29|) . The long-dashed curve is the 
result of Eq. (|3l)f derived following the ansatz @ . The 
short-dashed curve is our old result Q calculated follow- 
ing Eq. (fT9|) . We see that all three results (especially 
those calculated following Eqs. (j2U)) and ) are close 
to each other for temperatures T ^ 150 MeV. For higher 
temperatures deviations become more pronounced. 

In Fig. [21 the ratio of the bulk viscosity to the entropy 
density is presented as a function of temperature for two 
values of the baryon density ub = no and 4no, where 
no = 0.16 fm -3 is the nuclear density at the saturation 
point. The notation is the same as in Fig. [TJ The results 
are shown for the SHMC model (a) and for the ideal gas 
(IG) model with the same hadron set (b). We see that 
the curves calculated by Eq. (j2U]l and Eq. (|3"6"]l are closer 
to each other than to the curve calculated with Eq. (|19l) . 
For the IG, the curves calculated following Eqs. ([29)) 
and ([3"6"| coincide. Also, these curves are close to those 
estimated according to Eq. for T > 100 MeV. Com- 
paring figures (a) and (b) we also see that the presence 
of the quasiparticle interaction is more significant for low 
temperatures (T < 100 MeV) and it becomes less impor- 
tant for higher temperatures. 

Conclusion. Due to the complexity of the hadron 
system formed in actual heavy-ion collisions, it is hard 
to calculate the viscosity coefficients from the first prin- 
ciples with the help of the Kubo formulas. One could 
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FIG. 2: (Color online) The ratio of the bulk viscosity to the 
entropy density as a function of temperature for two values of 
the baryon density ns/no — 1 and 4 (from top to bottom). 



use the Kadanoff-Baym kinetic equations for the hadron 
resonances to derive general expressions for the kinetic 
coefficients, but at present realistic calculations do not 
seem possible even in the relaxation time approxima- 
tion 



12|. 



Therefore, making use of the quasiparticle 
Boltzmann-like equations, being treated within the relax- 
ation time approximation, can be considered as a forced 
step for practical evaluations of the kinetic coefficients in 
the given problem. The scaling hadron mass-couplings 
model of Ref. @ is an appropriate model for the de- 
scription of the equation of state of the hot and dense 



hadron system. The knowledge of the latter is necessary 
in order to perform evaluations of the kinetic coefficients. 
However, even an application of simplified phenomeno- 
logical expressions for the relaxation times of different 
species does not allow one to proceed in calculation of the 
bulk viscosity without doing additional assumptions, in 
particular the Landau-Lifshitz conditions should be ful- 
filled. However, these conditions cannot be satisfied on 
the class of solutions of the Boltzmann equations for our 
multi-component system treated within the relaxation 
time approximation, and additional ansatze are needed. 
Contrary, the result for the shear viscosity proves to be 
rather robust to these reductions, see @- 

We studied three ansatze previously used in the lit- 
erature and derived three expressions for the bulk vis- 
cosity (fl~9)). ([29)) and (|36|) . generalized to the case of 
nonzero chemical potentials ^b,s- Luckily, numerical 
evaluations performed following all three expressions de- 
viate not much from each other and confirm earlier re- 
sults • Although these results can be considered only 
as rough estimations, among them the result, Eq. (|29[) . 
seems to be the most theoretically justified. In order to 
perform more accurate calculations, one should go be- 
yond the scope of the relaxation time approximation and 
to fulfill the Landau-Lifshitz conditions on the class of 
solutions of the kinetic equation. However, such calcula- 
tions are much more involved than estimations presented 
in the given work and have not yet been carried out for 
multi-component systems with strong interactions. 
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